Conditional preparation of a non-classical state in the continuous variable regime: 

theoretical study 
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We study the characteristics of the quantum state of light produced by a conditional preparation 
protocol totally performed in the continuous variable regime. It relies on conditional measure- 
ments on quantum intensity correlated bright twin beams emitted by a non-degenerate OPO above 
threshold. Analytical expressions as well as computer simulations of the selected state properties 
and preparation efficiency are developed and show that a sub-Poissonian state can be produced 
by this technique. Projection onto a given trigger value is studied and then extended to a finite 
band. The continuous variable regime offers the unique possibility to improve dramatically the 
preparation efhciency by choosing multiple selection bands and thus to generate a great number of 
sub-Poissonian states in parallel. 

PACS numbers: 42.50 Dv, 42.65.Yj 



I. INTRODUCTION 

Conditional state preparation opens the possibility to 
generate a wide range of non-classical states of light. 
Quantum correlation between two modes -a signal and a 
trigger- is the prerequisite for such a general procedure. 
Pre-assigned events on the trigger condition the readout 
of the signal and cause a non-unitary state reduction. 

Various protocols have been proposed and imple- 
mented successfully in the photon counting regime. For 
instance, this strategy has been widely used to generate 
single photon Fock state. The required two-mode corre- 
lation can be obtained by an atomic cascade 1] or by the 
more efficient technique of parametric down conversion 
0. The single-photon state is created from the initial 
two-mode state when a photocount event is recorded in 
the trigger path. Two output ports of a beam splitter 
can also provide the required correlation: it has been 
shown theoretically that a squeezed vacuum transmitted 
through a beam-splitter is reduced to a Schrodinger-cat- 
like state when a photon is detected on the reflected port 
0. Non classical photon states have also been experi- 
mentally prepared by triggering on a given atomic state 

i- 

Recently, several conditional preparation protocols 
were introduced for the continuous variable regime where 
a continuously varying photocurrent is measured instead 
of single photon counts. For instance in continuous 
measurements are triggered by photocounts. The trig- 
gering condition and the characterization of the selected 
state can also be done both in the continuous variable 
regime. Some theoretical protocols have been proposed 
relying on two-mode squeezed vacuum generated by an 
optical parametric amplifier Q . 

The scheme we propose relies on quantum intensity 
correlated bright twin beams -signal and idler- gener- 
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FIG. f: Proposed conditional measurement protocol in the 
continuous variable regime. A ND-OPO generates above 
threshold orthogonally polarized and quantum intensity cor- 
related bright twin beams which are detected by high effi- 
ciency photodiodes. The continuously varying signal pho- 
tocurrent is kept only when the idler photocurrent takes a 
given value or falls inside a narrow band. 



ated by a non-degenerate optical parametric oscillator 
(ND-OPO) pumped above threshold 0. The procedure, 
sketched in Fig. 1, is the following : among all the 
recorded values of the signal beam intensities, one only 
keeps the events occurring in the time intervals when 
the idler intensity takes a given pre-assigned value, or 
more precisely, lies around this value within a small in- 
tensity range. All other time intervals are discarded. In 
the present paper, we study theoretically this conditional 
preparation technique. Its experimental implementation 
and the corresponding results have been detailed in 

This paper is organized as follows. In section^ we de- 
scribe the quantum state generated by a ND-OPO above 
threshold. Conditioning on a given value of the idler 
intensity results in a state reduction which generates a 
sub-Poissonian state. This procedure, presented in sec- 
tion lllll is then extended to the case of a selection interval 
of bandwidth A around a given value in section Hvl The 
state reduction is characterized according to the initial 
correlation, individual noise on each beam and position 
and width of the conditioning band. Section Ivl proposes 
a Monte-Carlo simulation of the reduced state properties. 
Section IVTI is then devoted to multiple-band selection, a 
possibility to significantly increase the efficiency of the 
conditional strategy which is only offered by schemes in 
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the continuous variable regime. In section rVIII the main 
conclusions of the paper are summarized and possible ex- 
tensions to more exotic non-classical states are discussed. 



II. QUANTUM STATE PRODUCED BY A 
ND-OPO ABOVE THRESHOLD 

The starting point of the theoretical analysis of the 
present conditional measurement is the precise knowl- 
edge of the quantum state produced by a non-degenerate 
OPO, from which the conditionally prepared state will be 
deduced by state reduction. When the twin photons are 
produced by spontaneous parametric down-conversion, 
the quantum state of the system is well-known and given 

by 



OO 

= (coshA)-i^(tanhA)"| 



(1) 



n=0 



where the indices 1 and 2 refer respectively to the sig- 
nal and idler modes and (resp. a^) is the annihila- 
tion (resp. creation) operator of a photon in mode i. A 
is proportional to the pump amplitude and crystal non- 
linearity. The state |ni,n2) is a Fock state with ni pho- 
tons in the signal mode and n2 photons in the idler mode. 
When the parametric down-conversion efficiency is very 
weak, which is experimentally the most frequent case, the 
state of the system reduces to the simple entangled state 



|*)~|0,0) + A|1,1) 



(2) 



When a single photon is detected in the idler mode, the 
system collapses by the state reduction process into the 
signal mode single photon state as is well known. 

To the best of our knowledge, there is so far no com- 
plete theory giving the state of light produced by an OPO 
in the Schrodinger picture, comparable to the Lamb the- 
ory giving the density matrix of the light produced by 
a laser The OPO has been described instead in 

most theoretical approaches in the Heisenberg picture, 
in terms of operators, which is not useful in the present 
case. Only Graham et al. have considered this prob- 
lem, mainly to determine the phase diffusion properties 
of the OPO, and not its joint photon-number distribu- 
tion. 

In the absence of a complete theory, we will use 
here a semi-phenomenological approach, relying on both 
theoretical considerations and experimental observations 
0, 13 • OPOs produce above threshold phase-coherent 
light in both signal and idler modes which have a Gaus- 
sian statistics and present close to threshold a large ex- 
cess noise compared to the shot noise level. We will call 
F the Fano factor of this photon distribution defined as 
the noise variance normalized to the one of a Poissonian 
distribution of same mean power: close to threshold, this 
factor is large. Furthermore, one knows that in the ab- 
sence of losses and in the case of a single output mirror 



of the OPO cavity, there is a perfect intensity correla- 
tion between the two modes when they are measured on 
time intervals long compared to the cavity storage time 
[l^. From this, we infer that the quantum state describ- 
ing the output of such a perfect OPO is the eigenvector 



of Ii — I2 ~ a^ai 



d2 with a zero eigenvalue. By 



assuming a gaussian photon statistics, the state can be 
described by 



|*OPo) = ^ Cn\n,n) 



n=0 



with 



1 



exp 



{n - nf 
2Fn 



or more generally by a density matrix 



n,n' 



(3) 



(4) 



(5) 



with Cn,n equal to the expression of |c„p given by equa- 
tion This form of the density matrix ensures at the 
same time a Gaussian distribution of photon numbers in 
each mode and a perfect intensity correlation between 
them. It will be shown later that the unknown expres- 
sions of Cn,n' for n ^ n' are not relevant in our calcula- 
tions. 

Uncorrelated losses in the system inevitably degrade 
the initial correlation between the signal and idler modes. 
These linear losses, that we will assume to be equal for 
the signal and idler modes, can be modelled by a beam- 
splitter of intensity reflection coefficient R equal to the 
loss coefficient and of transmission coefficient T = 1 — R 
(Fig. 12). If the input state impinging on such a beam- 
splitter is the tensor product \n, 0) of a Fock state on 
the input mode 1 and the vacuum in the loss mode, the 
output state l^'out) is 0| 



p=0 



with An^p given by 



An,p 



p\{n — p)\ 



T^~PRP 



(6) 



(7) 



corresponding to a binomial distribution of the photons 
over the two outputs of the beam-splitter, characteristic 
of a partition process. The system is then described by a 
pure state belonging to a four-mode Hilbert space : the 
signal, idler modes and the two loss modes of the two 
outputs. By tracing over these two loss modes which are 
not measured, one obtains the following reduced density 
matrix 

n' — \^ r , A A* A , A* 

H / ^ ^n.n' ■"-n,p^n' ,p-"-'n.,p' -"-n' ,p' 

n,n' ,p,p' 

\n — p,n — p')(n' — p,n' ~ p'\ (8) 
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Source of perfect 
intensity con'elated beams 
with gaussian statistics 
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FIG. 2: The dotted rectangle models a ND-OPO which can 
be described by a source of perfect intensity correlated beams 
plus a partition process generated by two symmetric beam- 
splitters with intensity transmission T and reflection R. The 
"gemellity" G of the beams is found to be equal to the loss 
coefficient R. 



This expression is valid for any values of the parameters 
i?, n and F. From it, one can calculate the probability 
distribution of the intensity difference. Its mean value is 
zero. Its variance can be written as 



= El^"l' E \An,r>\'\A^,p'\\p-p') 

n—0 p,p'—0 



(9) 
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and becomes 



A^{fi - 72) = 2RnT = 2Rn 



(10) 



where fiT ~ n' is the mean photon number of the signal 
and idler mode after the beam-splitter. A^(/i — I2) is 
independent of F and coincides with the value obtained 
by the usual linearized theory for the fluctuations. From 
this, one finds that the "gemellity" G, which is the re- 
maining noise on the intensity difference normalized to 
the total shot noise level 2n', is equal to the loss coeffi- 
cient R. 

As bright beams have very large mean values in com- 
parison with the Poissonian standard deviation, one can 
easily assume that both the reflected and transmitted 
beams have also large mean values, even for R ot T close 
to 1. Thus, one can use the Stirling's approximation and 
show that, in a normalized form 



|A 



1 



\/27rni?T 



exp 



(p - nRf 



2nRT 



(11) 



Let us calculate the photon distribution of the individ- 
ual beams after the beam-splitters. If we only consider 
measurements involving output 1, the reduced density 
matrix p" is obtained by tracing over the output 2 

^ c„,„'A„,pA;, p|n-p)(n'-p| (12) 

n,n' ,p 

The normalized probability distribution for the number 
of photons can thus be written as 



P(ni) = (ni|p"|ni) 



(13) 



and becomes 



+ CXD 



Pini) = E 



-i \/27rFn yj2TmRT 



(14) 



xexp 



(riT - n'f (n-ni- nRf 



2Fn'T 



2nRT 



We introduce 5 — nT — n' and e = ni — n' and we 
assume that n can be taken constant and equal to n in 

the factor . As the numbers of photon are large, 

the parameters 6 and e can be considered as continuous 
variables. The sum can thus be rewritten as an integral 
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2R{S + rV) 
Developed to the first order in Jr and 4j 

1 1 



dS 



P{s) 



exp 



1 



T y/2TiFn y/2Trn'R 
1 1 



(16) 



' 2Fn'T 

and by using the formula 

f +00 



2Rn' ' Rn' 



2Rn' 



d6 



/ exp [-AiS'^ + A2S] dS ^ . — exp 
J -00 V Ai 



A. 



(17) 



the probability distribution can be integrated to give 
1 r [ni 



P{ni 



V2nn'F' 



exp 



2n'F' 



with 



F' ^ R + FT 



(18) 



(19) 



n=0 



Thus we find that the new Fano factor F' after the beam- 
splitter is given by i? + FT. In the case of important 
losses (R ^ 1, T ^ 0), this Fano factor goes to 1, i.e. 
the output state distribution is Poissonian as expected. 



III. STATE REDUCTION BY A PHOTON 
NUMBER MEASUREMENT 

In this section, we study the reduced state resulting 
from a photon number measurement on beam 2 giving 
the pre-assigned value N. By tracing over the trigger 
mode, the reduced density matrix p'" for the signal mode 
only can be written as 

\n-p)(n'-p\ (20) 
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One can see that the photon probabihty distribution de- 
pends only on the diagonal terms of the density matrix 
p. The photon probability distribution of beam 1 when 
a given value N is measured on beam 2 is thus 



+ OC 



I ^n.n- 



N\ 



(21) 



n=0 



With F' and n' the Fano factor and mean value estab- 
lished in the last section, this distribution can be read 
as 



P{ni,N) 



1 



1 



V2^ ^ 2-KnRT 

^ n— 



exp 



{nT - 



n 



l\2 



(nT~nif {nT~Nf 



2nRT 



2nRT 



2{F' - R)n' 
(22) 



Let us consider e — ni — n' introduced in the previous 
section and a = N — n' which corresponds to the distance 
between the value of conditioning N and the distribution 
center n' . Using calculation techniques similar to the 
ones of the previous section, the probability distribution 
can finally be written as 



P(e,a) 



1 



^/2^Tn'F' 
1 



exp 



a 



exp 



2n'F' 
2n'Vr 



with 



Vc = 2G 



G_ 
F' 

F' 



(23) 



(24) 



(25) 



Vc is the conditional variance of the intensity fluctua- 
tions the signal beam knowing the intensity fluctuations 
of the trigger beam. This parameter plays an important 
role in the characterization of Quantum Non Demolition 
measurements .14] . 

The first factor between parenthesis in equation H23() 
expresses the preparation probability, i.e. the proportion 
of selected points in an experimental implementation. It 
is maximum when a = 0. A value chosen in the wings 
of the gaussian distribution leads to a smaller efficiency. 
The preparation probability is inversely proportional to 
V F'n' . This factor is thus very small in the case of bright 
beams: for 1064 nm beams with a mean power of 1 mW 
and a Fano factor equal to 100, the maximal preparation 
probability obtained when a = is as small as 5.10"^". 

The second factor of this probability shows that the se- 
lected state has a gaussian photon distribution centered 
around the value /?a, slightly different from the trigger- 
ing value a except for a = 0. This difference can be 
interpreted by the dissymmetry of the probability distri- 
bution around the conditioning value. For a conditioning 
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FIG. 3: Conditional variance as a function of the Fano factor 
for different initial intensity correlation between signal and 
idler beams, measured by the "gemellity" G. Values of G as 
small as 0.1 have been obtained experimentally [l^ . 



value equal to the mean of the distribution (a = 0), this 
difference vanishes as the distribution is symmetric. 

According to H23|l , the reduced state has a photon vari- 
ance equal to the conditional variance Vc- This shows 
that the selection process we have used in the conditional 
measurement has extracted the maximum information 
available from the correlation of the two beams and has 
transferred it to the signal beam. This beam can exhibit 
a sub-Poissonian photon distribution when Vc < 1. Fig- 
ure shows the conditional variance as a function of 
the Fano factor for different gemellities G. If the corre- 
lation is perfect, this protocol generates a number state. 
It is worth noting that for a gemellity inferior to 0.5 the 
distribution is sub-Poissonian whatever the initial Fano 
factor. For large Fano factors, which is the experimental 
case, the noise reduction is twice the initial gemellity as 
Vc can be approximated by 2G. 

The conditional variance Vc also characterizes the noise 
reduction obtained when an active feed-forward correc- 
tion of signal beam intensity is implemented by opto- 
electronics devices controlled by information collected on 
the idler The two techniques which take advantage 
of the quantum correlation have therefore the same ulti- 
mate performance. 



IV. STATE REDUCTION BY A BAND 
SELECTION 



As shown in the previous section, triggering on a sin- 
gle value of photon number in our regime of bright beams 
where f\! is very large leads to a close to zero preparation 
probability incompatible with experimental implementa- 
tion. An interesting question is to determine the photon 
distribution when one projects onto a finite band instead 
of a given value. It is worth pointing out that it is always 
the experimental case because of the limited acquisition 
precision. 

Let us consider a conditioning band of width A around 
the value N = n' -\- a. The photon probability distribu- 
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FIG. 4: Photon probability distributions of conditionally pre- 
pared state. The probability distributions are all normalized 
to 1 at their maximums for clarity of the plots. The horizon- 
tal unit is the width a = \ffr! of the Poisson distribution of 
same mean intensity. Thick curves give the initial gaussian 
distribution and dashed ones correspond to a coherent state 
of same mean intensity. Arrows give the extension of the se- 
lection band, (a) and (b): the selection band is centered on 
the mean value (a = 0). The bandwidth is A = O.lcr for 
(a) and lOo" for (b). A large band results in a wide distri- 
bution (Kurtosis coefficient lower than 3). (c) and (d): the 
band is centered on a = lOcr. A = O.Ict for (c) and IOct for 
(d). Non central selection band results in an asymmetry of 
the distribution (non null skewness) (F' = 100, G = 0.18) 



tion of the reduced state can now be written 



noise distribution of the initial state and of a coherent 
state of same mean intensity have been superimposed. 
One observes a narrowing of the probabihty distribution 
below the shot noise level in the case of a very narrow 
bandwidth for any band center. When the bandwidth 
increases, the reduced state differs more and more from 
a gaussian distribution. 

The distance from a gaussian distribution can be char- 
acterized by the coefficients of skewness and kurtosis 
which are related to higher moments of the probability 
distribution The coefficient of skewness provides 

a dimensionless measure of the asymmetry of the prob- 
ability distribution. Skewness is null when the band is 
centered on the mean value. For a normalized probability 
total surface, the kurtosis coefficient distinguishes distri- 
bution which are tall and thin from those that are short 
and wide. It has the value of 3 for a gaussian distribu- 
tion and a value lower than 3 corresponds to a less peaked 
distribution. Let us note that the third order cumulant 
corresponds to the kurtosis excess. A very small value 
of kurtosis is associated to a square-shaped distribution. 
A perfect gaussian distribution is obtained in the limit 
of very narrow bandwidth -as calculated in the previous 
section- or in the case of very large bandwidth which 
corresponds to keep all the values and the selected state 
corresponds thus to the initial state. A narrow band of 
selection results in a close to gaussian distribution with 
a kurtosis coefficient almost constant around 3. This co- 
efficient decreases when the band is extended (Fig. O. 
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and integrated as 
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FIG. 5: Kurtosis coefficient of the conditionally selected state 
according to the selection bandwidth normalized to ct = 
(b): zoom for small bandwidth, (a = 0, F' = 100, G = 0.18) 

We can derive from equation l|27|l an approximate but 
analytical expression for the probability distribution by 
expanding in powers of 



where erf denotes the error function. 

From this exact expression, one can plot the proba- 
bility distribution of the conditionally prepared state in 
different cases. The initial Fano Factor and gemellit 
are taken from our ND-OPO experiment detailed in 
[F' = 100, G = 0.18). Figure ^ gives the photon prob- 
ability distribution in different possible configurations. 
The band can be centered around the mean {a = 0) or 
around an arbitrary value taken in the wings of the ini- 
tial gaussian distribution (a ^ 0) . The bandwidth A can 
also be small or large relative to the standard deviation 
cr of a Poisson distribution of same mean intensity. The 
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24n'T/2 



(-1 + 




with P{e, a) the probability distribution determined in 
the last section. 

Let us underline that for A = 1 the photon distribution 
is equal to the one established in the last section. This 
value corresponds to the minimum bandwidth necessary 
to discriminate two consecutive photon-number states. 



FIG. 6: (a) Intensity noise on the reduced state and (b) 
preparation probability (proportion of selected points) as a 
function of the selection bandwidth normalized to a. (a = 
0,F' = 100, G = 0.18) 



The third order term in A can be neglected in this 



approximation as long as A ^ a 
probability can then be written 



n' . The preparation 



N. 



selected 



total 



(AM 
V2^' 



exp 



{a/al 
2F' 



(29) 



where Ntotai and Ngeiected correspond respectively to the 
total number of events and to the number of selected ones 
by the conditional protocol. 

As a result, conditioning on a narrow band results in 
a preparation probability proportional to A without de- 
creasing the non classical character of the projected state. 
This range of A is the optimal one in order to generate 
a sub-Poissonian state. With a bandwidth equal to O.lcr 
the efficiency reaches 0.4%. 

From equation (|27|l . it is possible to compute the ex- 
act properties of the state for any value of A. In Fig. 
(0, we give the noise reduction and the probability of 
preparation for a narrow bandwidth. One can see that 
for very narrow bandwidth the preparation efficiency in- 
creases linearly with the bandwidth whereas the squeez- 
ing is almost constant. This corresponds to the first order 
in A as appears in H28I) . For a larger band, the efficiency 
still increases but at the expense of a decreased noise 
reduction. The efficiency would be higher for an initial 
state with less excess noise above shot noise. This case 
can correspond to a ND-OPO pumped well above thresh- 
old where it has been shown theoretically that both twin 
beams can be individually squeezed 18] and still highly 
correlated. In the case of narrow bandwidth, the band 
center has no effect on the non classical character but re- 
sults in a lowered preparation efficiency due to the gaus- 
sian distribution of the initial noise (Fig. (TJ). It is 
worth noting that the large extension of the initial noise 
should permit to implement a great number of indepen- 
dent selection bands. Section lVll is devoted to this partic- 
ular conditional strategy offered by continuous variable 
regime. 
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FIG. 7: Preparation probability as a function of the center 
of the narrow A = O.lcr selection band. {F' = 100, G = 0.18) 



sion. However, our protocol is easy to test by Monte- 
Carlo simulations. 

In order to implement the protocol, one needs to pre- 
pare two random arrays A and B with super-Poissonian 
statistics and exhibiting a given amount of correlations. 
These two arrays will correspond to the fluctuation dis- 
tribution of the signal and idler twin beams. Actually, 
we generate three independent arrays. Let us call C a 
super-Poissonian distribution and X and Y two Poisso- 
nian distributions. If the correlation between the signal 
and idler beams is perfect, we associate to the signal and 
idler beams the same distribution C . For a finite amount 
of correlations, a similar strategy to the one detailed in 
section is performed. We add to the initial common 
distribution C independent contributions which can be 
seen as vacuum contributions when the initial state is 
incident on a beam splitter with reflection coefficient R. 
The statistics distribution of signal and idler after the 
partition process can be defined by 



B = y/RY + Vl- RC 



(30) 



The parameter R is the same as the one used in previous 
sections and determines the gemellity of the beams. 

The selection of relevant events can then be done. It is 
worth noting that conditional measurement experiments 
are very frequently made after the end of the physical 
measurement, i.e. by post-selection of the events Q- So, 
when the two previous distributions have been generated, 
we dispose of the same kind of data that after an exper- 
iment. 

We have simulated all the properties presented before. 
As expected, the simulations are in perfect agreement 
with our theoretical predictions. Such simulations also 
should allow to test the protocol when n' is not too large, 
a case where the approximations leading to equation H23() 
are no longer valid. 



V. MONTE-CARLO SIMULATION OF THE 
STATE REDUCTION 



VI. MULTI-BAND SELECTION 



So far we have examined the conditionally prepared 
state from its theoretical probability distribution expres- 



In any conditional protocol, the preparation probabil- 
ity in a given time interval is an important parameter 
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FIG. 8: Experimental multi-band selection. The horizontal 
unit is the width a of the Poisson distribution of same mean 
intensity. Independent bands of width A = 0.2a have been 
implemented. Band are separated by 4a in (a) and 2a in (b). 
{F' = 100, G = 0.18) 

to characterize its efficiency. The efficiency of all con- 
ditional preparation techniques is usually very low. Im- 
proving the brightness of the source is thus an important 
and topical challenge in the photon counting regime . 
The preparation probability is also low in our protocol 
as shown in Fig. © and (0 but the continuous detec- 
tion results in a great number of selected points in a very 
short recording time Q . 

Furthermore, the continuous variable regime offers the 
possibility to improve the efficiency of the conditional 
measurement strategy by a large factor. This possibil- 
ity does not exist for the single photon counting case. 
One can implement multiple selection bands with differ- 
ent centers on the idler intensity. Independent selection 
bands will correspond to independent sets of time win- 
dows. In each of them, the state is reduced to a given sub- 
Poissonian state with a constant noise reduction, wher- 
ever the band center. By using independent intervals, 
one keeps most of the values of the idler intensity and so 
improve by a large factor the success rate of the prepa- 
ration. This scheme -directly related to the nature of 
continuous variable- opens the possibility to generate in 
parallel a great number of sub-Poissonian states. Figure 
© gives the experimental implementation of multi-band 
selection for different distances between the band center. 
Experimental details are given in Q. 

VII. CONCLUSIONS 

In the photon counting regime, conditional measure- 
ments plays a crucial role as illustrated by the various 
schemes theoretically proposed or experimentally imple- 



mented. The extension to the continuous variable regime 
is of great importance as the efficiency can be dramati- 
cally improved and more complex triggering schemes im- 
plemented. 

The required quantum correlation is generated in our 
protocol by a ND-OPO above threshold. It has been 
shown that conditioning on a given value of the idler in- 
tensity leads to the preparation of a sub-Poissonian state 
with a noise reduction equal to the conditional variance, 
i.e. the initial gemellity minus 3 dB in the case of a 
large Fano factor. We have extended this triggering con- 
dition to a finite band. For a band largely narrower than 
the standard deviation of a coherent state of same mean 
intensity, the preparation probability increases linearly 
with the bandwidth and the noise reduction remains al- 
most constant. As opposed to the discrete variables case, 
there is a trade-off between the preparation efficiency and 
the non classical character of the selected state. Theo- 
retical calculations and computer simulations are in very 
good agreement, and both account very well for the ex- 
perimental results detailed in 8]. 

The next step should be to extend this conditional 
measurement strategy to the preparation of more ex- 
otic non-classical states as it is the case in the photon 
counting regime. One can think for instance to gener- 
ate Schrodinger-cat-like state. It could be also of great 
interest to extend this protocol to EPR beams, in par- 
ticular generated by a self-phase-locked OPO where the 
anti-correlated phase fluctuations of the twin beams are 
accessible poj . 
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